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Abstract 

We systematically investigate general forms of mass matrices for three-generation up 
and down quarks, including asymmetrical ones in generation space. Viable zero matrix 
elements are explored which are compatible with the current observation of masses and 
mixing angles, and also with the recent measurement of CP violation in the B-meson 
system. The simplest form with the maximal number of vanishing matrix elements is 
found to be almost consistent with the experimental data, but has a scratch that one 
of the mixing angle is slightly large. At the next-to-minimal level, it is found with a 
help of leptonic generation mixing that only six patterns of mass matrices well describe 
the experimental data. These sets of mass textures predict all the properties of quarks, 
including the CP violation, as well as the large (charged) lepton mixing, which may be 
appropriate for the atmospheric neutrino in grand unification scheme. 



1 Introduction 



It is certain that one of the most important issues confronting the standard model is the 
generation structure of quarks and leptons. Various neutrino oscillation experiments have 
recently been bringing out the generation structure of the leptonic sector. The Super- 
Kamiokande experiment has established the neutrino oscillation in the atmospheric neu- 
trinos with nearly maximal mixture pQ. As for the solar neutrino problem, the Mikheyev- 
Smirnov-Wolfenstein solution 75} is strongly suggested by the recent experimental results |2j 
if there exists the neutrino flavor mixing between the first and second generations 4 . On 
the other hand, the other mixing (the 1-3 mixing) has been found to be rather small 
similarly to the quark sector. While these experimental progresses have been giving us a 
new perspective beyond the standard model, it seems that we are far from a fundamental 
understanding of the origin of fermion masses and mixing angles. 

A promising approach to the issue of the generations is to assume that some of Yukawa 
matrix elements are vanishing. An immediate and important consequence of this approach 
is to reduce the number of free parameters in the theory and to lead to relations among the 
fermion masses and mixing angles. Moreover that would provide a clue to find symmetry 
principles or dynamical mechanisms behind the Yukawa sectors, which are interpreted 
as remnants of fundamental theory in high-energy regime. Most of the previous work 
along this direction, including the systematic analysis by Ramond, Roberts and Ross [5], 
assumed that the matrices of Yukawa couplings are symmetric about the generation indices 
(hermitian matrices).* However, in the context of the standard model and even in grand 
unified theory, it is not necessarily required that Yukawa matrices are symmetric. In fact, 
the present experimental data indicate that the (large) leptonic mixing mentioned above 
is quite different from the (small) quark mixing. It is interesting that this asymmetrical 
observation is known to be compatible with quark-lepton grand unification if the fermion 
mass matrices take asymmetrical forms in generation space (Sj. There are also some classes 
of non- hermitian ansatze for quark mass matrices jS], which are consistent with the data 
and cannot be transformed to the solutions obtained in Ref. It is therefore worthwhile 
to do systematical examination and to complete the classification of viable asymmetric 
forms of fermion mass matrices. 

In this paper, we investigate phenomenologically viable mass matrices of up and down 
quarks, assuming that the Yukawa couplings generally take asymmetric forms in generation 
space. In particular, we look for as simple forms as possible, that is, mass matrices with 
the maximal number of vanishing elements. Vanishing matrix elements are expected to be 
deeply connected with underlying physics, such as flavor symmetries, in more fundamental 
theory to shed some lights on constructing realistic models of quarks and leptons. Note 
that our treatment is general and includes symmetric mass matrices as limited cases. 

This paper is organized as follows. In Section 2, we describe the Yukawa sectors of 
up and down quarks in the standard model and introduce the parameterization needed in 

* Systematic studies and classifications of fermion mass matrices and mixing have also been performed 
in other approaches \7\. 
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later discussion. Our analysis does not depend on any details of the Higgs field profile, and 
therefore can be straightforwardly applied to other cases such as grand unified theory and 
supersymmetric models. Sections 3 and 4 are devoted to analyzing which forms of matrices 
(vanishing matrix elements) are compatible with the current experimental data. It is found 
that the minimal (Section 3) and next-to-minimal (Section 4) cases contain only a few 
types of mass matrices phenomenologically viable. We summarize our results in Section 5. 

2 Formulation 

In this section, we briefly review the formulation of the up and down quark Yukawa sectors 
in the standard model. The SU(S) x SU(2) x U(l) gauge-invariant Yukawa interactions 
are given by 

-C Y = Q l (Y u ) lj u Rj H* + QiiYjijdRjH + h.c, (2.1) 

where Qi denote the SU(2) doublets of left-handed quarks, and ur, cIr are the right-handed 
up- and down- type quarks, respectively. The Yukawa couplings Y u and Yd are 3x3 matrices 
the generation indices), and H is the SU(2) doublet Higgs field. After the electroweak 
symmetry breaking, these Yukawa interactions lead to the following quark mass terms: 

-C m = u L i{M u )ijU R j + d L i{M d )ijd R j + h.c, (2.2) 
(M u )ij = (Y u ) ijV , {M d ) i:j = (Y d ) ijV , 

where v is a vacuum expectation value of the neutral component of the Higgs field H. If 
the model is super symmetrized, the Yukawa terms are described by superpotential in terms 
of quark and Higgs superfields. Only one difference is that in the supersymmetric case two 
types of Higgses must be introduced to have gauge-invariant Yukawa terms. This procedure 
does not bring any modification to the structure of Yukawa couplings Y Ujd , and therefore 
the following analysis is straightforwardly extended to supersymmetric models and also to 
other scenarios. 

The generation mixing is physically described by the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix which consists of two unitary matrices 

Vcku = vl L V dL . (2.3) 

These unitary matrices diagonalize the mass matrices M u and M d ; 

M u = V uL M°V^ R , M d = V dL M°Vj R . (2.4) 

The diagonal elements of and correspond to the experimentally observed mass 
eigenvalues. 

With phase degrees of freedom of the six quark fields, the number of observable pa- 
rameters in the CKM matrix is reduced to four (the overall phase rotation is physically 
irrelevant). However, it is important to distinguish the contributions of V u l and V d L from 
a viewpoint of pursuing clues to find more fundamental theory of quarks and lepton such 
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as grand unification and flavor symmetry. A generic 3x3 unitary matrix U has 9 free 
parameters and can be parameterized as 

U = $Oi$'0 2 03*"- (2-5) 
The matrices O, (i = 1, 2, 3) represent the rotations in the index space around the i-ih axis 



1 1 \ / cos 6 2 sin0 2 \ / cos# 3 sin0 3 0^ 

Ox= cos Ox sin 6»i , 2 = 1 , 3 = -sin0 3 cos0 3 
\0 -sin6»i cos 0i/ \-sin0 2 cos0 2 / \ 1, 

(2-6) 

The diagonal phase matrices 3>'s are given by $ = diag.(e J9i , e 11 ^, 1), <J>' = diag.(l, 1, e tuJ ) 
and $" = diag.(e jp , e iq , e lv ). There generally exists three rotation angles and six complex 
phases. When applied to the above quark mixing matrices V u l and V^l, the phases factors 
in $" are always unphysical degrees of freedom since they can be absorbed by field redefi- 
nitions of the quark mass eigenstates. It will also be found that phase matrices $"s do not 
appear throughout this work (except for only a few examples discussed at the beginning 
of Section 4) because we will consider 3x3 matrices with non-vanishing determinants and 
at most five independent elements. In this CclSG, Si matrix M is always expressed such 
that M = JM T J' where J and J' are the diagonal phase matrices and M r contains only 
real parameters. Thus the matrices MM^ and M^M are diagonalized by real orthogonal 
matrices, up to overall phase rotations corresponding to <!> or in (|2.5|) . On the other 
hand, the phases factors in $ of the up and down sectors generally contribute to the CKM 
matrix elements. Thus the CKM mixing matrix is found to be written as 

^ckm = OlPO d , (2.7) 
P = K$d, Oi = O u 2i Oz h 

where the subscripts i = u,d label the up- and down-type quarks, respectively. As men- 
tioned above, two complex phases in P play an important role for reproducing CP-violating 
quantities since their changes generically affect the CKM matrix elements. The numerical 
results of such phase factors will be discussed in later sections. 

The experimentally observable quantities in the quark Yukawa sector are 3 mixing 
angles with 1 complex phase in the CKM matrix and 6 mass eigenvalues. The measured 
values of the CKM matrix elements are jlUj 



|Vckm| 




'0.9739 - 0.9751 0.221 - 0.227 0.0029 - 0.0045\ 
0.221 - 0.227 0.9730 - 0.9744 0.039 - 0.044 
0.0048 - 0.014 0.037 - 0.043 0.9990 - 0.9992 / 

(2.8) 

The various experimental observations of CP-violating phenomena yield the CP violation 
in the standard model, which are translated to ^U] 

J C p = (2.88 ±0.33) x 10~ 5 , (2.9) 
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where Jqp denotes the reparameterization- invariant measure of CP violation This 
value corresponds to the Kobayashi-Maskawa phase in the standard parameterization as 
<5km = 60° ± 14°. Moreover the recent results of studying the decay of the B mesons to 
charmoniums indicate j!2j 

sin 20i //3 = 0.726 ±0.037, (2.10) 

where (f>± = j3 is one of the angle of the unitary triangle for the 5-meson system, which 
is defined as (f>i = /3 = arg(V* d V c b/V^ d V t b)- This is the angle that is most precisely known 
at present and is expected to provide the most stringent constraint. The current-quark 
masses at the ^-boson mass scale are evaluated |131 including various effects such as the 
QCD strong coupling factors and we obtain 

m u = 0.000975 - 0.00260, m d = 0.00260 - 0.00520, 

m c = 0.598 -0.702, m s = 0.0520 -0.0845 (2.11) 

m t = 170 - 180, m b = 2.83 - 3.04, 

in GeV unit. The most recent results of the Tevatron CDF and D0 experiments indicate 
the top quark mass which is a bit larger than that quoted above |14j . However, if taken 
into account, the analysis of mass matrix forms presented in this paper is not significantly 
changed, since the most influential ingredients are practically the masses of lighter gener- 
ations. In the following analysis, we use these experimental data as input parameters and 
explore possible forms of quark mass matrices. 

3 The Minimal Asymmetric Matrices 

We would like to systematically search for the mass matrices of up and down quarks which 
are consistent with the current experimental data. Our analysis is based on possible zero 
elements in the mass matrices. Namely, the aim of this paper is to investigate how small 
number of non-vanishing matrix elements can account for the existing data. Here the 
number of zeros means independently-vanishing elements in a mass matrix. In particular, 
for a symmetric matrix, "1 zero" implies that a diagonal element or a pair of off-diagonal 
elements in symmetric positions takes a negligibly small value. 

Let us first consider the mass matrix of up- type quarks. In the present work we assume 
that the up-quark mass matrix M u is symmetric. This assumption is motivated by grand 
unified theory, where the left- and right-handed up quarks in one generation often belong 
to the same multiplet of unified gauge symmetry, like SU(5) and larger. In this case one 
obtains the identical mixing matrix for left- and right-handed up quarks; V u l = V U R = V u . 

It is first noticed that three independent, non- vanishing matrix elements are needed 
to reproduce the observed mass eigenvalues of the three-generation quarks. Moreover a 
determinant of mass matrix must be nonzero. The minimal forms of matrices which satisfy 
these criterions are found to coexist with at most three zeros, and the independent matrices 
are given by the following three types; 
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• Mul 

M u = | b | (3.1) 

• Mu2 

M u = | 6 | (3.2) 

• Mu3 

M u = la bj (3.3) 

The matrix elements a, b, c are nonzero and the blanks denote vanishing entries. All other 
forms of matrices consistent to the criterions can be obtained by relabeling the generation 
indices. For example, by exchanging the first and second generations (both for the left- 
and right-handed fermions), the matrix Mu2 is converted to the form discussed in |15j 

M u = a 
\b 

It should be noted that we do not assume any hierarchical orders among the non-vanishing 
elements. Therefore the analysis of the above three types of matrices (Mul, Mu2, and 
Mu3) includes the whole possibility of symmetric matrix with three zeros. The matrix M u 
is diagonalized as 

(771 \ 
-m c Vl (3.4) 

m t ) 

The negative sign in front of m c is just a convention introduced in order that the mass 
eigenvalues and the parameter c are real and positive. With a suitable phase redefinition 
of the up- type quarks, we take the non- vanishing matrix elements to be real parameters 
without loss of generality. These values can be fixed by the three mass eigenvalues from 
the following three equations: 

tr M u = m u — m c + mt, (3-5) 
trM^ = ml + m 2 c + mf, (3.6) 
det M u = —m u m c m t . (3-7) 

Thus a unitary matrix V u which diagonalizes a three-zero symmetric M u is described in 
terms of the up-quark mass eigenvalues. 

For the down sector, the mass matrix is not necessarily symmetric. The minimal crite- 
rion for a realistic mass matrix is the same as for symmetric matrices; three independent, 
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non- vanishing matrix elements, and a non- vanishing determinant. To satisfy these require- 
ments, we have at most six zero elements. A matrix with six zeros has only three parameters 
which correspond to three mass eigenvalues. That only gives no mixing angle or exchang- 
ing generation indices. However since any type of up-quark mass matrices with three zeros 
cannot be diagonalized with the observed CKM matrix, a matrix with asymmetric six zeros 
is not suitable for the down sector. We thus find that the most economical candidates for 
a realistic mass matrix of down-type quarks have five zeros. They can generically describe 
three eigenvalues and one mixing angle. It is found that there are 36 types of mass matrices 
with five zeros and non-vanishing determinants. At this stage, since we are not requiring 
any hierarchy among matrix elements as in the case of up-quark mass matrices, these 36 
(=6x6) patterns are related to each other through the permutations of three rows and/or 
three columns. Namely, one can obtain all the patterns by exchanging the generation labels 
from a single matrix, e.g., 



The matrix elements d, ■ ■ ■ ,g are made real- valued by phase redefinitions of quark fields. In 
our convention, a permutation of columns corresponds to a rotation of generation indices of 
the right-handed quarks, which rotations do not change the mass spectrum and the CKM 
matrix elements. On the other hand, the exchanges of rows, i.e. relabeling three left-handed 
down quarks, do affect on the observable mixing angles. This is because we have already 
used the label exchange degrees of freedom to reduce the number of matrix patterns for the 
up-type quarks. Therefore all possible permutations of rows must be taken into account in 
the down sector to explore the whole combinations of up- and down-quark mass matrices. 
We thus consider the following 6 types of mass matrices as the minimal candidates with 
five zeros; 




(3.8) 



• Mdl 




(3.9) 



• Md2 




(3.10) 



• Md3 




(3.11) 



• Md4 




(3.12) 
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• Md5 



Ma 



'd 

e 

k9 fj 



Md6 



M d 




(3.13) 



(3.14) 



As expected, these 6 patterns are transformed to each other by changing the generation 
indices of the left-handed quarks, up to permutations of the right-handed ones. 

Given the possible forms of mass matrices, Mul-Mu3 and Mdl-Md6, we analyze which 
combinations of mass matrices explain the observed masses and mixing angles. It is easily 
found that the matrices Mul and Mu2 cannot fit the data. This is because they have 
zero or one finite mixing in the up sector and all the candidates of Md can induce only 
one generation mixing, that necessarily results in the CKM matrix with more than one 
vanishing entries. The only remaining possibility is the matrix Mu3 for the up sector. The 
matrix which diagonalizes Mu3 is approximately written by the mass eigenvalues 



/ l / flu 



\Ou\ 



m c j m u \ 



J 



1 



(3.15) 



This shows that the 1-2 mixing angle from the up sector is roughly given by ^Jm u /m c = 
0.037 — 0.066. Therefore in order to generate the observed Cabibbo angle, a 1-2 mixing 
angle from the down sector is required to be of order O(10 _1 ), which selects out Md3 
or Md4 for an appropriate matrix for the down quarks. The matrices Md3 and Md4 are 
diagonalized by rotations of the first and second generations and do not affect the third 
column of Vckm- Consequently the combinations (Mu3, Md3) and (Mu3, Md4) predict the 
CKM angles involving the third generation 



\V ub \ 
\Vtd\ 

\Vcb\ 




m t 



V, 



as 



0.00233 - 0.00391, 

0.00894 - 0.0122, 



TO,, 



m. 




Vt,\ ^ J — 




m c 

0.0576 - 0.0643. 



m t 



(3.16) 
(3.17) 
(3.18) 



It is found that the first two predictions well agree with the experimental values (|2.8j) but 
the mixing between the second and third generations is slightly larger than the observation. 

We have found in this section that the simplest and realistic forms of quark mass 
matrices can accommodate symmetric three zeros and asymmetric five zeros in the up and 
down sectors, respectively. Their explicit forms are highly constrained by the experimental 
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values of mass eigenvalues and mixing angles, and there exists only two possibilities which 
are given by 

<d g 

Ma = | a h | . My, - ( e ) . (.119) 

and 

M u = [a b | , M d = \ g e | . (3.20) 

ft 

It is clear that a simultaneously exchange of the identical generation indices for ul, ur, and 
di completely preserves the physical consequences. Moreover, as noted in the classification, 
any permutation of the right-handed down quarks (i.e. of the columns of Md) is also allowed 
phenomenologically. Noticing this fact, one can see, for instance, that the matrix M d in 
(|3.20j) reconciles the Georgi-Jarlskog ansatz ^f)], up to unphysical field rotations, and could 
easily be extended to including the charged-lepton mass matrix. A numerical evaluation 
for the combination (|3.19|) presents us an example 

0.000221 (A 5 - 56 ) \ 

M u = ( 0.000221 (A 5 - 56 ) 0.05 78 (A 1 - 88 ) J m t , (3.21) 

0.0578 (A 1 - 88 ) 0.997 (A - 00 ) / 

0.00518 (A 3 48 ) \ 

M d = j -0.0285 (A 2 - 35 ) \m b , (3.22) 

1.00 (A - 00 )/ 

where the couplings are chosen so that we can fit as many observables as possible to the 
experimental data. A non-trivial phase factor is also required, e.g. P = diag.(e -0 ' 7?n , 1, 1), 
in order to reproduce the observed value of CP violation (|2,9|) . In this case, however, we 
find that alternative indication of CP violation (|2. 1U|) cannot be reproduced. The unitary 
triangle for the l?-meson system is distorted due to a large value of \V c b\, while the area 
of the triangle is correct. A similar result is obtained for the combination (|3.2()|) since 
physical consequences are now determined modulo right-handed mixing of down quarks. 
In the above example, we have alternatively written down in the parentheses the exponents 
of a small parameter A (= 0.22) for non-vanishing elements in M u and M d . Such expressions 
in terms of an expansion parameter would be suitable to gain an insight into the fermion 
masses problem in a view of flavor symmetries. 

Though there is one unsatisfied point for the above two combinations of mass matri- 
ces, i.e. a slightly large value of the 2-3 CKM mixing angle {\V c b\ — y/m c /mt), one can 
easily find some remedies. The discrepancy may be removed with radiative corrections, 
for instance, the renormalization-group effects between the electroweak and high-energy 
scales. A probable source of such effects is a flexibility of the top-quark Yukawa coupling 
in high-energy regime due to its fixed-point behavior in the infrared. To see how the 
renormalization-group evolution makes the situation better, let us first define the following 
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ratio: 

The problem with the mass matrices (j3.19|) and (|3.2U|) is that their prediction A pre ~ 1 is not 
consistent to the experimental results A ex p > 1.31 at the electroweak scale. Within a good 
approximation that the third- generation Yukawa couplings are dominant, the evolution of 
the mass ratio m c /mt and the mixing angle V c b is governed by the equations 

din (2*) -1 r 

""A = i6^[C(^33) 2 + ^33) 2 ], (3-24) 

^ - ^[(W + C*»fl, (3.25) 

where t = In [i denotes the renormalization scale. The direct contribution from gauge 
couplings is generally irrelevant to the running of mass ratios and mixing angles. The 
coefficients £ and rj are model-dependent constants, for example, £ = 3/2 and rj = —3/2 in 
the standard model, and £ = 3 and rj = 1 in supersymmetric standard models (though the 
1^33 effect is negligible in the standard model). Thus we obtain 

^ = 3^2 ~ 0(W + v(Y d33 ) 2 ] . (3.26) 

It is noted that the coefficient in front of (Y u33 ) 2 is negative in usual Higgs doublet models. 
This negative sign suggests that the ratio A in high-energy regime is reduced from the 
value observed at the electroweak scale. Therefore the renormalization-group evolution in 
fact ameliorates the problem with the mass textures (j3.19|) and (|3.20|) . It is also noticed 
that, since the second term in the right-handed side of (|3.26|) is positive for supersymmetric 
cases, a smaller value of ^33 is preferred to cure the mismatch between A pre and A exp . In 
the case that only the top Yukawa coupling is dominant, we obtain by integrating (|3,26|) 
over the range between /i = 1 TeV and the unification scale A, 

A(A) _ fY u33 (A)\^fg 1 (A)\^fg 2 (A)\^fg 3 (A)\^ ^ QQ / F w33 (A) \^ 



\Y u33 (fj,)J \gi(fJ,)J XSbitiJ V 53(^)7 V y «33(^). 
for the minimal supersymmetric standard model. In the second equation, we have roughly 
assumed that the three gauge couplings of the standard gauge groups 51,2,3 are unified at 
A in one-loop order estimation. The situation is more improved for the standard model 
because of a larger negative value of the coefficient 2r\ — £ and the negligible positive 
contribution from the bottom Yukawa coupling, and one typically obtains 



-1 



A (A) /Y u33 (A)\^ 
K J - 0.59 " 3J ; J . (3.28) 



a 00 ' VWm). 

at a high-energy scale A ~ 10 14 GeV. It is clearly seen that \V c b\ and ^Jm c jmt at high- 
energy regime become closer than around the electroweak scale, and the texture ansatz 
(|3.19j) and (j3,2()|) would work better in high-energy theory such as grand unified models. 

On the other hand, a more direct resolution to the problem A pre 7^ A exp is to incorporate 
additional non- vanishing elements into the Yukawa matrices. This is the option we will 
explore in the next section. 



9 



4 The Next-to-Minimal Asymmetric Matrices 



The previous analysis has shown that, at classical level, any combinations of symmetric 
three-zeros M u and asymmetric five-zeros Md are too simplified to be totally consistent 
with the observed data. In this section we investigate a possibility to relax the constraints 
on matrix forms and to introduce one more non- vanishing matrix element. 

The first case to consider is to work with symmetric two zeros in the mass matrices of 
up quarks. An interesting observation is that, when the 2-2 element in M u is turned on, 
the five-zeros matrices Md in 1)3.19)) and 1)3.20)) may completely explain the data. This is 
because the mixing angle V c b could be controlled by a free parameter in M u , irrespectively 
of the charm quark mass, which resolves the difficulty discussed in the previous section. 
The zero structures of such mass matrices are ruled out at the 3a level by the current 
experimental data if non- vanishing elements are symmetric and hierarchical valued )17j . 
Such a special case is obtained from our general form by exchanging the first two indices 
of right-handed down quarks in 1)3. 2U|) and identifying the 1-2 and 2-1 matrix elements. 
Up to relabeling generation indices, there are 4 types of mass textures with symmetric two 
zeros. Exploring all possible patterns for M u (with symmetric two zeros) and Md (with 
asymmetric five zeros), we find that the following mass matrices are successful to explain 
the experimental data: 



M d 



a 
d 
b 



h 
f 



b 
d 



f 



9/ 



All 4 = 2x2 combinations of M u and Md well describe the present experimental data. 
It should be noted that there are additional solutions with relabeling the generation in- 
dices, while physical implications are unchanged. First, any permutation of dn (i.e. of the 
columns of Md) is phenomenologically allowed. In addition, simultaneously exchanges of 
the identical generation indices for ul, ur, and cLl completely preserves the physical con- 
sequences and also become the solutions. The above sets of mass matrices are consistent 
with all the properties of quarks, including the recent measurements of CP violation in the 
£>-meson system ()2.10)l . 

The other case we will pursue in the following is to extend the down-type mass matrix 
Md to contain (asymmetric) four zeros. It is found that there are 81 types of matrices 
with nonzero determinants, which are related through the permutations of rows and/or 
columns. All possible forms can be generated from the following four representatives of 
mass matrices: 



'd 



'd 



'd 



e 
9 



d 



(4.1) 
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It is first noticed that the last matrix generates only one generation mixing angle and does 
not cure the problem in the previous section, where one of the CKM mixing angles is not 
entirely consistent to the experimental data. Therefore we can safely drop this matrix (and 
the other 8 ones generated by changing the labels) in the analysis below. For the former 
two mass matrices in (|4.1|) . a permutation of the first two columns produces the same 
modification as that obtained by permuting the first two rows, since any hierarchical order 
among the matrix elements is not supposed at this stage. This fact reduces by half the 
number of independent forms of matrices generated by exchanging the generation labels 
from the first two representatives in (|4.1jl . The total number of independent mass matrices 
we will consider is 72 (= 18 + 18 + 36). 

As in the analysis of Section 3, if one does not count up the rotations of the right-handed 
down- type quarks (the columns of matrices), only the following 12 types of matrices should 
be taken into account: 




(4.2) 



Since we now drop the degrees of freedom of the right-handed down rotations, physical 
consequences can be read from the symmetric matrix M^mJ, which has 6 independent 
elements. For the above 12 types, a hermitian matrix M^mJ contains 5 free parameters and 
necessarily leads to one vanishing element or one relation among non-vanishing elements. 
Thus an asymmetric four-zeros matrix gives similar results to that of a symmetric one- 
zero matrix, as far as mass eigenvalues and left-handed down mixing are concerned. It is 
noticed that there is a difference between M^mJ and symmetric mass matrices discussed 
in jS] that none of the diagonal elements of M^mJ can be zero for matrices with non- 
vanishing determinants. It was found in the analysis of Ref. |B] that symmetric two-zeros 
mass matrices for the down sector are consistent with the experimental data. Given these 
facts, the above asymmetric four-zeros mass matrices are expected to account for the proper 
mass eigenvalues and mixing, because they correspond to symmetric one-zero matrices 
which have one more free parameter. The number of combinations which can explain the 
data is hence rather large, and the exploration along this line unfortunately seems not to 
provide a new perspective for the origin of fermion masses and mixing angles. 

Let us proceed the discussion taking account of the roles played by the mixing of right- 
handed down quarks. In the standard model, the mixing of right-handed (5C/(2)-singlet) 
fermions is irrelevant to the CKM mixing and unphysical (unobservable) degrees of freedom. 
This is not necessarily true in various extensions of the standard model. For example, in 
supersymmetric extensions of the standard model, the right-handed mixing of fermions 
which diagonalizes a Yukawa matrix is transfered to that of corresponding scalars via 
supersymmetry-breaking scalar masses. Thus the masses of scalar superpartners generally 
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have generation dependences and cause observable effects, such as flavor-changing decays of 
heavy fermions. A more interesting situation arises in the frameworks of grand unification. 
In this case, quarks and leptons are unified into some multiplets of unified gauge group 
and the mass matrices of quarks are often closely related to those of leptons. This fact 
may give rise to an apparent difficulty in simultaneously realizing the small CKM mixing 
and the observed large lepton mixing, which is described by the Maki-Nakagawa-Sakata 
(MNS) matrix Vmns [TH] - The parallelism between quarks and leptons, which is a sign 
of grand unification, does not seem to work in the Yukawa sector. There is however an 
interesting observation that the mixing angles of left-handed charged leptons are correlated 
to those of right-handed down quarks and therefore the CKM mixing does not necessarily 
connect with the MNS mixing. This idea is easily achieved in SU(5) grand unification and 
larger unified theories [SJ, where a key ingredient is that an anti 5-plet of SU(5) contains 
one- generation right-handed down quark and left-handed lepton doublet, and also there 
are some multiplicities of anti 5-plets, which are naturally incorporated in SO(10) or Eq 
unified models. This mechanism automatically makes asymmetric since it changes the 
property of anti 5-plets only, while keeping that of 10-plets. 

In this way the mixing of right-handed fermions is not necessarily unobservable quan- 
tities. In the following analysis, motivated by grand unification view mentioned above, we 
investigate possible connections between V^r and the leptonic mixing matrix which diag- 
onalizes the mass matrix of three- generation charged leptons. In particular, we examine 
which combinations of mass matrices for the up and down sectors suggest large leptonic 
mixing, recently observed in various neutrino experiments. Since the flavor rotation of 
right-handed down quarks is now supposed to change physical consequences, it is not used 
to reduce the number of candidates for four-zero matrices as done in the above. We will 
therefore exhaust all of the most generic 72 candidates for realistic down-quark mass ma- 
trices with four vanishing entries. 

Assuming that the neutrino oscillations account for the solar and atmospheric neu- 
trino data, the recent experimental results indicate rather large angles for the 1-2 and 
2-3 generation mixing in the MNS matrix, but a small one for the 1-3 mixing angle 
04lNs)i3 < 0.14 — 0.22 L 5_. As for the two large mixing angles, the best fit value for 
the atmospheric neutrino data is the maximal mixing (9 ~ vr/4), and on the other hand, 
the solar neutrino deficit needs a large but non-maximal value of the 1-2 mixing angle [I]. 
In the following analysis, we first consider, just as a first approximation, the maximal angles 
both for the 1-2 and 2-3 mixing, and then examine possible deviations from these maximal 
angles. The lepton mixing matrix is defined as 

Vmns = V^V V , (4.3) 

where V e L rotates the three-generation charged leptons such that the mass matrix of charged 
leptons is diagonalized, and V v denotes some mixing matrix in the neutrino sector. Its form 
crucially depends on the neutrino property and we leave it, together with detailed analysis 
of neutrino mass texture zeros, to another future task [T2J|. It should be noted that this does 
not mean that we take V u = 1 in the following analysis. In fact, our result will show that 
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considerable contribution to lepton mixing is needed to come from neutrino mass matrices, 
which could be realized in a huge variety of neutrino models. 

As we mentioned, if embedding the theory into grand unification scheme, the mixing of 
charged leptons V b l may be related to that of right-handed down quarks as V e L — VdR, up to 
corrections due to the breaking of unified gauge symmetry. To precisely reproduce the mass 
eigenvalues of charged leptons, it is in fact needed to take in some breaking effects which 
split the properties of quarks and leptons. Typical examples of such splitting are provided 
by the Georgi-Jarlskog factor ^fj] and higher-dimensional operators involving Higgs fields 
that break quark-lepton symmetry. We assume, just for simplicity, that such breaking 
effects are small not to significantly change the analysis below. We are thus interested in 
the following typical forms of mixing matrices for the down sector, which are associated 
with large generation mixing of left-handed charged leptons:^ 

no o \ 

V dR = 1/V2 -1/V2 , (4.4) 

\o 1/V2 1/V2 J 

'l/y/2 -1/V2 0\ 
V dR = ( 1/V2 1/V2 j , (4.5) 

-1/2 1/2 \ 

V d R = I 1/V2 1/2 -1/2 . (4.6) 

\ 1/V2 1/V2J 

A general procedure for examining viable forms of VdR is as follows. At first, evaluate 
the matrix 

M\M d = V dR | ml V} R , (4.7) 

where the matrix VdR is parameterized as given in Section 2. We consider in this section 
the matrices Md with asymmetric four zeros. Namely, they contain five free parameters, 
and a matrix of the form M^M has six independent elements. The equation 1)4. 7 jl therefore 
imposes one constraint which can be used to eliminate a mixing angle of the right-handed 
down quarks. As noted above, one or two additional constraints are obtained to reduce 
the number of independent (mixing) parameters, when one explores the solutions of VdR 
with large mixing. For example, we can fix (V^h)32 = l/v2 m the case of 1)4.4)1 . Once the 
matrix elements in Md are solved with respect to remaining independent parameters, the 
mixing matrix for left-handed down quarks is expressed as 

V dL = M d VdR ( mj 1 J . (4.8) 



^The mixture l|4.5|l would be excluded by the neutrino oscillation experiments as it would generate too 
a large value of the 1-3 lepton mixing angle, if the atmospheric neutrino angle comes from the neutrino 
sector. While included into the analysis, as we will show, the mixing pattern 114.511 is already disfavored by 
the quark data alone. 
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Such left-handed down mixing is used to examine which forms of mass matrices produce 
the observed values of the CKM matrix elements. 

Numerically exhausting all the possible forms of mass matrices, we find that the down- 
quark mass matrices with asymmetric four zeros are unfavorable to a sizable value of 
mixing angle between the first and second generations of right-handed down quarks. This 
fact is easily understood in the case that there is large mixing only between the first two 
generations. Such a mixing matrix is given by 

/cos6> -sin0 0\ 
V dR = sin/9 cos<9 (4.9) 
V l) 

at the leading order of other small mixing angles. The case where 9 ~ 7r/4 corresponds to 
Eq. ([4.5)1 and could provide a solution to the solar neutrino problem in grand unification 
schemes. To explain the quark mixing angles for the up-quark matrices Mul, Mu2 and 
Mu3, the left-handed mixing in the down sector needs to satisfy 

/ 1 0(e) 0(e 3 )\ 
VdL ^ 0(e) 1 0{e 2 ) , (4.10) 
\0(e 3 ) 0(e 2 ) 1 / 

where e is a small parameter of order of 10 _1 . It is found from the analysis in the previous 
section that the marginal requirements are sizable contributions from the down sector to 
V us and V& (and not necessarily to V u b). This is translated to lower bounds on the left- 
handed mixing of down quarks, for example, | (V^^)i2 1 > 0.16 and | (V^^)23 1 > 0.012. When 
there is a solution for the above-described procedure, the corresponding down-quark mass 
matrix is given by 

'm d 

m s I I - sin 9 cos6> I . (4.11) 





It is clearly seen that four zeros in Md cannot be realized since the matrix elements in 
the second and third rows are always non-vanishing for any precise values of VdL (|4.10|) 
satisfying the lower bounds mentioned in the above. The situation might be improved by 
turning on fluctuations around the exact form of VdR (|4.9|) . In this case, one is in fact able 
to take either {Md)zi or {Md)z2 as zero, if additional 0(e 4 ) mixing in V d R is introduced. 
However not all of the elements in the first row become zero, in particular, either (Md)u 
or (Md)i2 can be set to zero. We thus find that the down-quark mass matrices with 
asymmetric four zeros generically conflict with large 1-2 mixing in the right-handed down 
sector. 

It turns out that solutions with two large mixing like (|4.6j) are also absent. In the limit 
of bi-maximal mixing of d,R, the down-quark mass matrix becomes 

'm d 

M d I 0(e) 1 0(e 2 ) II m s J ^ § -4= | . (4.12) 
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The situation is rather different from the case of (|4.11|) . for instance, it is now possible to 
have vanishing matrix elements in the second row in (|4.12j) due to the presence of the third- 
generation large mixing. Note that the third-row elements are necessarily non-vanishing, 
even if one introduces sizable deviations from the maximal or zero mixing angles in V d R, 
i.e. / 7r/4 and/or / 0. Accordingly it is enough to consider physical consequences 
of the matrices M d with four zeros placed in the first and second rows (and the other 
components are nonzero) [2U]. If one adopts the up-quark mass matrices Mul-Mu3, the 
mixing angles from the down sector have some lower bounds, in particular, | (Vrf^)i2 1 > 0.16 
is needed. It is numerically evaluated that the condition |(Vrf^)i2| > 0.16 constrains the 
other mixing angles as \(V dR ) 2 i\ < 0.40 (0.44) for \(V dR ) 32 \ = 0.7 (0.61). This value is 
translated, in the limit of negligible 1-3 mixing, to the upper bound of the 1-2 mixing 
angle #3 < 26.1°, which is excluded at more than 3er level by the recent results of neutrino 
experiments, if generation mixing in the neutrino sector is found to be small. We thus 
find that any four-zero mass matrix in the down sector is not compatible with bi-large 
generation mixing of right-handed down quarks. 

Finally let us consider the case that the right-handed down mixing between the second 
and third generations is large and the others are suppressed [Eq. (|4.4j) ]. This implies, 
if adopting the grand unification, the amount of mixture of the atmospheric neutrinos. 
Following the general procedure described before, we have exhausted the possible patterns 
and found that, at classical level, the following mass matrices satisfy the criterion for 
(charged) lepton mixing, while the quark masses and the CKM matrix elements (including 
the KM phase) are properly reproduced: 



M d 



9 

e 



All 6 = 2x3 combinations of M u and M d are consistent to the present experimental data. 
Note that one type of mass texture for M d [the second one in (|4.2|) ] almost describes the 
data we have listed in Section 2. It is however found that the whole parameter space of 
that texture is excluded by the measured value of CP violation in the iT-meson system (the 
€k constraint). All the other combinations of mass textures are not compatible to the ex- 
perimentally allowed parameter region for the atmospheric neutrino problem at more than 
6<7 level, unless there is sizable contribution to mixing angles from the neutrino sector. The 
examples of numerical fits for these matrix elements are shown in the appendix. It should 
be noted that there are additional solutions with relabeling the generation indices, while 
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physical implications are unchanged. The combinations obtained by exchanging the second 
and third columns of Md are viable. This is simply because we now consider the situation 
that the second and third generations of cIr are largely mixed. In addition, simultaneously 
exchanges of the identical generation indices for ul (and ur for symmetric textures) and 
di completely preserves the physical consequences and also become the solutions. No other 
exchanging symmetry does not exist. The above sets of mass matrices are consistent with 
all the properties of quarks, including the recent measurements of CP violation in the B- 
meson system, as well as large lepton mixing for the atmospheric neutrino problem. More 
conservatively, they provide sizable contribution to leptonic 2-3 mixing, while satisfying 
the experimental results of quark masses and the CKM matrix elements. It is interesting 
to see that, in the down sector, the numerical exploration of parameter space shows that 
a correlation f — g should hold for all the above solutions. The above list of textures 
contains the Fritzsch ansatz 21,, but M u is symmetric while is not [221) leading to 
large generation mixing in VdR- As mentioned in the beginning of this section, such an 
asymmetric form of mass matrix often plays a key role for neutrino physics in grand unified 
theory We have shown that the above three forms of are the minimal extensions 
of the ansatz / ~ g to including the first generation. The obtained matrices are successful 
to explain the observed quark masses and mixing angles and have the maximal number of 
vanishing elements. 

5 Summary and Discussions 

The study of the origin of fermion masses and mixing angles is one of the most important 
unresolved issues in particle physics. As a plausible approach to this issue, possible zero 
elements in mass matrices have been extensively examined and the obtained results have 
suggested useful guides for realistic model construction. In this paper we have systemati- 
cally investigated what types of quark mass matrices with non-symmetrical forms can be 
consistent with the experimentally obtained CKM matrix and mass eigenvalues. Our first 
principle is that a mass matrix has as simple form as possible, namely, to search for the 
minimal number of free parameters in the mass matrices. This leads us to consider some of 
mass matrix elements to be vanishing. The existence and structures of zero matrix elements 
are expected to be deeply connected with underlying physics, such as flavor symmetries, 
in more fundamental theory of quarks and leptons. We have first examined experimentally 
viable mass matrices in the case where the up-quark sector has symmetric three zeros and 
the down-quark sector asymmetric five zeros. This is the simplest possibility apparently not 
to conflict with the experimental data, and can almost explain the observed quark masses 
and mixing angles. The situation is rather different from the case where the down-quark 
mass matrix contains at most four vanishing elements. We then find that there exist var- 
ious forms of mass matrices consistent with the existing experimental data, and it seems 
difficult to find some clues to understand the generation structure. Additional informa- 
tion comes from the recent observation of neutrino generation mixing. If working with the 
grand unification hypothesis, the mixing of S'C/(2)-doublet leptons is correlated to that of 
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SJ7(2)-singlet down-type quarks. To investigate the implications of large mixing angles in 
the lepton sector, we have searched viable solutions which induce large right-handed mixing 
in the down sector, and found that there only exist six patterns of mass matrices with a 
large mixture between the second and third generation of right-handed down quarks. Fur- 
thermore it turns out in our framework that the large angle solution for the solar neutrino 
problem cannot be realized from the charged lepton sector with asymmetric four- zeros Md 



The observed large amount of mixing angle of solar neutrinos then should come from 
the neutrino sector. If the minimality principle is applied to neutrino mass matrices, 
the simplest matrix forms (i.e. with the maximal number of zero matrix elements) could 
be found out. However the neutrinos have rich phenomenology and their property has 
not been fixed experimentally. In particular, there still exists wide possibilities for the 
neutrino mass spectrum, which fact generally makes the thorough analysis of neutrino mass 
textures laborious. We here briefly discuss several results for possible forms of neutrino mass 
matrices which lead to the large lepton mixing between the first and second generations 
and have the maximal number of allowed zero matrix elements. First, consider the effective 
neutrino mass operator KijLiLjH* H where Lj denote the three-family lepton doublets. 
The coefficient matrix Kij is symmetric in the generation space. This higher-dimensional 
operator induces the Majorana neutrino mass matrix Ml = k(H*H) after the electroweak 
gauge symmetry breaking. We find two types of the simplest forms of Ml which have 
symmetric four zeros and are given by 



It is interesting to note that this form of neutrino mass matrix predicts the spectrum with 
the inverted mass hierarchy and an exactly massless neutrino for the third generation. 
Moreover, taken into account the observed neutrino mass differences, such Ml leads to 
almost maximal mixing angle between the first two generations. As a result, it can be 
matched with only four of six combinations of quark mass matrices found in the previous 
section. If the minimality analysis is extended to the next-to-minimal level, i.e. Ml with 
symmetric three zeros, we find thirteen patterns are allowed, each of which predicts char- 
acteristic mass spectrum of light Majorana neutrinos. Another well-known scheme for light 
neutrinos is to consider (3 x 3) Dirac neutrino mass matrix M v and right-handed Majorana 
one Mr. Also in this case, the analysis is quite different from the quark sector, mainly 
because one neutrino can be massless. We find from the exhaustive exploration that the 
maximal number of vanishing matrix elements is ten which consists of asymmetric seven 
(six) zeros in M v and symmetric three (four) zeros in Mr. There are four patterns of the 
seven-zeros M u cases, which contain as an example 



(M e ). 




(5.1) 




(5.2) 
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and six patterns for the six-zeros M v cases, for example, 



M v = 




Mr 



V 



v 



w 



) 



(5.3) 



All the ten patterns of M v and Mr generate light neutrino mass matrix Ml in the form 
of ()5.1|) after the seesaw operation. Therefore the resultant mass spectrum and possi- 
ble partners for quark mass matrices are the same as the cases 1)5. 1J) . Some different 
phenomenology may appear through lepton flavor-violating processes induced by lepton 
Yukawa couplings The detailed analysis of minimal lepton mass matrices and their 

phenomenological implications will be presented in a separate paper |19| . 

In the analysis of this paper, except for the discussion at the end of Section 3, we have 
not taken into account of the dependence of matrix forms on the renormalization scale, and 
considered generic features of 3x3 quark Yukawa couplings including asymmetrical matri- 
ces. For more precise treatment, the renormalization-group evolution of Yukawa couplings 
are required to be evaluated, because zeros of matrix elements should be implemented at 
some high-energy scale such as a grand unification scale. The observable quantities at the 
electroweak scale deviate to some extent from the values estimated in high-energy regime. 
However one of the most important points is that the fermion mass ratios of the first to 
second generations is almost insensitive to radiative corrections due to the fact that the 
dominant contribution to flavor-changing evolution comes from the Yukawa couplings of 
the third generation. In our analysis, the selection of viable forms of mass matrices has 
mainly depended on whether the down-quark matrices satisfy the experimental value of 
the 1-2 CKM mixing in conspiracy with the up sector. It is therefore expected that the 
renormalization-group analysis does not destabilize the results of our analysis of possible 
zero elements, while there certainly exist some scale dependences of non-vanishing matrix 
elements in the presence of significant contributions from the gauge and top-quark Yukawa 
couplings. This latter fact is supposed to only change 'initial' values of non-vanishing 
Yukawa couplings at a high-energy scale. The results presented in this paper are also use- 
ful for explaining the flavor structures of quarks and leptons in grand unification schemes. 

Finally, we would like to comment on some phenomenology related to the solutions 
obtained in Section 4. These solutions predict similar sizes of off-diagonal elements to 
the 3-3 elements and radiative corrections from Yukawa couplings are important for flavor 
physics. For example, if the theory is supersymmetrized, flavor violation in the Yukawa 
sectors is translated to off-diagonal components of supersymmetry-breaking scalar masses 
through the radiative corrections. That could induce sizable rates of flavor-changing neutral 
currents for quarks and charged leptons in supergravity models. Further searches of 
flavor-violating processes will provide us a new perspective of flavor structures in high- 
energy regime. 

In the viewpoint of distinguishing possible solutions, it is important to examine ob- 
servable signals of underlying theory. In addition to signals of underlying symmetries or 
dynamics, the improved measurements of low-energy observable quantities allow us to dis- 
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criminate discrete ambiguities of possible matrices. As for the solutions 1-6 presented in the 
appendix, it can be seen from the numerical analysis that the solutions 3 and 4 have sizable 
contributions to (VdR) 31 components. This means that they predict (Vmns)i3 ~ O(l0~ 1 ) if 
there appears no fine tuning of parameters in V^r and V u . Since the planned improvements 
in the sensitivity to (Vmns)i3 are expected to reach 0.05 (23], these solutions would be 
supported or disfavored when a precise value of (Vmns)i3 is measured. For other genera- 
tion mixing, the solutions 3-6 are found to have relatively larger values of (VdR) 21 (of the 
order of the Cabibbo angle) than the other solutions. This fact also distinguishes possible 
textures, for example, if the theory is extended to incorporate supersymmetry (breaking) 
or grand unification. Together with these issues stated above and others, it is hoped to 
find what underlying theory governs the masses and mixing angles of quarks and leptons. 
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A The order estimation and numerical evaluations of the 
quark mass matrices 



In this appendix, we would like to present the order estimation of quark mass matrix el- 
ements and typical examples of numerical fitting for the solutions obtained in Section 4, 
where the up-quark matrices have symmetric three zeros and the down-quark ones asym- 
metric four zeros. Since the observed values of masses and mixing angles are hierarchi- 
cal, one could parameterize matrix elements by integer exponents of a small parameter A 
(= 0.22) times O(l) coefficients, which originate from the ambiguities of Yukawa coupling 
constants. Such expressions with integer exponents might be useful for getting ideas of con- 
structing fermion mass matrix models with flavor symmetries. We have found in Section 4 
that there are 6 = 2x3 combinations of up and down quark mass textures well describe 
the current experimental data. The order estimation of these mass matrix elements are 
presented in Tabled where we have not explicitly written down O(l) coefficients mentioned 
above. 



M u /m t 



A 8 



A 2 

A 2 A 



A 6 

A 6 A 2 
A 2 A 



A* v 



A 2 

A 2 A , 



A 6 

A 6 A 2 

A 2 A°; 



A' y 



A 2 

A 2 A 



A 6 

A 6 A 2 
A 2 A 



M d /m b 



f \ 4 









Table 1: The typical orders of matrix elements for six possible mass textures. We have 
not explicitly included O(l) coefficients, which would be needed to precisely reproduce the 
experimental data. Note that there are also additional solutions obtained (i) by exchanging 
the second and third generations of da (columns in M^) and/or (ii) by identically relabeling 
generation indices for ul, ur and d^. 



Suitably choosing the 0(1) coefficients (i.e. Yukawa couplings) in the textures listed in 
Tabled we obtain numerical examples for these six solutions (Table |2J). 
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1 


2 


3 


4 


5 


6 


m u 


0.00179 


0.00104 


0.00179 


0.00223 


0.00179 


0.000983 


m d 


0.00387 


0.00470 


0.00283 


0.00321 


0.00511 


0.00271 


m s 


0.0562 


0.0683 


0.0637 


0.0805 


0.0566 


0.0542 


m c 


0.613 


0.611 


0.652 


0.633 


0.601 


0.621 


m b 


2.90 


3.01 


2.94 


2.99 


2.91 


2.87 


m t 


178 


177 


175 


176 


179 


171 


\v us \ 


0.225 


0.225 


0.223 


0.223 


0.224 


0.223 


\V cb \ 


0.0433 


0.0430 


0.0433 


0.0433 


0.0428 


0.0434 


\v ub \ 


0.00403 


0.00414 


0.00426 


0.00429 


0.00377 


0.00437 


Jcp /10- 5 


2.90 


2.60 


3.16 


2.79 


2.91 


2.83 


sin 2^1 /p 


0.709 


0.692 


0.762 


0.735 


0.694 


0.742 


\(V dR )2l\ 


0.0116 


0.00990 


0.144 


0.167 


0.285 


0.135 


\(v dR ) 32 \ 


0.681 


0.682 


0.643 


0.642 


0.656 


0.707 


\(V dR )3l\ 


0.0108 


0.00922 


0.125 


0.147 


0.0136 


0.00918 



Table 2: Numerical examples for the predictions of the texture combinations given in 
Tabled The mass eigenvalues are denoted in GeV unit. 
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